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Abstract 

The high-peaked BL Lac object Pks2155-304 shows high variability at multiwavelengths, i.e. 
from optical up to Te V energies. A giant fl are of around 1 hour at X-ray and TeV energies 



was observed in 2006 lAharonian et al.l (120071) . In this context, it is essential to understand the 
physical processes in terms of the primary spectrum and the radiation emitted, since high-energy 
emission can arise in both leptonic and hadronic processes. In this contribution, we investigate 
the possibility of neutrino production in photo-hadronic interactions. In particular, we predict a 
direct correlation between optical and TeV energies at sufficiently high optical radiation fields. 
We show that in the blazar Pks2155-304, the optical emission in the low-state is sufficient to lead 
to photo-hadronic interactions and therefore to the production of high-energy photons. 
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1. Very high-energy emission from blazars 

The very first detection of a blazar at TeV photon energies h appened in 1992, when Mkn421 



was observed with the HEGRA telescopes dPunch et all 119921) . Since then, large progress has 



been made detecting more and more active galactic nuclei (AGN) with detailed spectral fea- 
tures. As of today, 26 blazars as well as 2 radio galaxies of Fanaroff-Riley typ^3 are known to 
emit at very high-energy (VHE, i.e. E y > 100 GeV). Blazars typically show a spectrum with 
a double-hump structure, i.e. a low-energy peak at radio to optical wavelengths (or X-ray) and 
a high-energy peak at X-ray (or gamma-ray) to TeV energies. The entire emission is believed 
to arise in the relativistic s hock fronts in AGN jets, where charged particles are accelerated to 
ultra high-energies (UHE) dMeli et allEooilBiermann et alll2009t e.g .). Then, the low-energy 



component is due to synchrotron radiation by the accelerated electrons. Several processes can 
produce VHE radiation: With relativistic electrons accelerated in the shock fronts of the blazar's 
jet, synchrotron photons are produced. The latter are scattered to VHE via the Inverse Compton 
process with the primary electrons. Hadronic Cosmic Rays (CRs) are believed to be accelerated 



1 The radio galaxies are M 87 and Cen A, see Wagner 1 2009) for a complete list of known sources 
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along with the electrons, up to 10 21 eV. These CRs can radiate at VHE through synchrotron radi- 
ation, or interactions with ambient matter or photon fields. Here, we consider the possibility of 
CR interactions with photon fields. 

Blazars are observed to be highly variable at all wavelengths. The source of t his va r iabilit y 
is still under discussion, rangin g from in stabilitie s in the acc retion disk, see e.g. IWiital d2006h . 
to precession of the jet, see e.g. lGergely and Biermann (2009). Other explanations are jets with 
density inhomogeneities or turbulent jets, see Romero (1995) and references therein. Such phe- 
nomena lead to the variation of the flux at low energies (between radio and X-ray). The variability 
at the higher energies (X-ray to TeV energies) depends on the VHE radiation origin: 

1. Inverse Compton scattering can arise in two different scenarios: On the one hand, it can 
be due to the interaction between the accelerated electrons with the synchrotron photon 
field ("Synchrotron Self Compton", SSC). This leads to a direct correlation between the 
low- and high-energy humps of the spectrum. On the other hand, if external photons are 
boosted to VHE by electrons, such a correlation is not necessary. The high-energy hump 
would then show a variability uncorrelated to the low-energy spectrum. 

2. Hadronic interactions can lead to the production of VHE photons through the production 
of 7T°, as neutral pions decay into two VHE photons. Here, a correlation between the low- 
and high-energy hump is not necessary. It can, however, be present if the hadronic CRs 
interact with the synchrotron photon field itself. In the case of proton-proton interactions, 
no correlation would be present. In addition, synchrotron proton radiation can contribute 
at VHE. 

3. A combination of hadronic and leptonic VHE emission is likely to be present in many cases, 
as protons and electrons are co-accelerated in the jet's shock fronts. The continuous inter- 
action of the different radiation fields and charged particle populations may lead to a series 
of loops of VHE production. This may lead to a strongly non-linear correlation between 
the intensity of different wave l ength s. A description of repeated leptonic loops is given in 



Kellermann and Paulinv-Toth ( 19691) . The same mechanism may work for hadronic loops 



or a combination of both, always assuming a low-energy target photon field. 

As the synchrotron field from primary electrons is a good target for Inverse Compton scattering, 
a direct correlation between low- and high-energy photons is usually interpreted as such. In 
so-called orphan flares, a strong variability at VHE is observed in the absence of low-energy 
variability. This in turn is often taken to be an indication as a result of hadronic interactions, as 
these do not require an intensity change at low energies. In this paper, we present a model where 
variability at low-energies leads to a correlated signal at high-energies. 



2. Photo-hadronic interactions in blazars 



When CRs are accelerated to UHE, they can interact with the present photon fields via the 
A-resonance: 
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The decay of the A-resonance leads to the production of charged an d neutral pions, which decay 
into high-energy neutrinos and photons, respectively (iBecken. 120081 e.g.)): 



yy. 
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The probability for n + production is 2/3, while n° are produced in 1/3 of the cases. Around 
1 /5 of the proton energy goes into the pion, and the decay products carry equal energy, so that a 
fraction of 0.25 (0.5) is transferred to each neutrino (photon). Furthermore, we need to consider 
that the neutrinos oscillate from 2 muon neutrinos and one electron neutrino to one electron, 
one muon and one tau neutrino. Considering these factors, the high-energy photon and neutrino 
spectra are connected a^| 
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Photo-hadronic interactions occur under the following conditions: 



1. Protons need to be accelerated to sufficiently high energies. 

2. The photon field must be dense enough for the protons to interact, i.e. we need an optical 
depth for proton-photon interactions of t P7 ~ 1 or larger. 

The first conditio n is believ e d to b e fulfilled, as AGN jets enable proton acceleration up to 
10 21 eV, see e.g. iMeli et al.l (12008b : iBiermann et al.1 (120091) . An appro ximation of the opti- 
cal de pth for photo-hadronic interactions in knots of blazars is given by iBecker and Biermannl 
(2009), 
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Here, L y is the luminosity observed from the blazar at a given frequency v. Further, F w 10 
is the Lorentz factor of the shock and 6 » 0.1 is the jet's opening angle. The parameter e^ ot 
gives the fraction of luminosity coming from a single knot in the jet. The value eknot = 0.1 indi- 
cates that the luminosity comes from 10 knots, which is the order of magnitude to be expected. 
Finally, Zj is the distance of the knot along the jet. The fi rst shock is expected to be present at 



around 3000 gravitational radii, r g , from the foot of the jet dMarscher et al.L 120081) . As the optical 



depth drops with the distance from the center and we need those strong shocks to have efficient 
neutrino production, this first strong shock provides the best environment for the production of 
secondaries. The flaring time scale of Pks2155-304 also suggests that the emission happens in a 
compact region, underlining the fact that it could be due to photohadronic interactions. 



3. Multiwavelength observations of Pks2155-304 

Pks2155-304 is a blazar at a redshift of z = 0.116, located in the southern hemisphere, at 
(RA, Dec) = (21h58'52.0 ", -30°13'32") in J200 coordinates dFev et all |2004 . It was the 5th 
source reported at y-rays (Cha dwick et al. , 1999), with a steady-state emission of around 20% of 
the Crab nebula. Pks2155- 304 shows a very rapid variability, with extreme flares, the most in- 
tense one observed in 2006 (lAharonian et all 120071) . with an increased intensity of more than one 
order of magnitude at TeV energies. This extraordinary flare lasted about 60 minutes and showed 
a peak intensity of 7 times the Crab nebula. CANGAROO-III w as triggered by the H.E .S.S. ob- 
servations and detected the flaring state at complementary times ISakamoto et al.l (12008b . During 
this flaring state, no other wavelength bands were observed. 



2 for a detailed description, see HalzerJ 120081) . 
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In a later flare, a multiwavelength campaign of Chand ra and H.E.S.S. was ar ranged to study 
the correlations between the different wavelength bands dCostamante et al. , 2008 ). Here, X-ray 
and TeV variabilities were more moderate, i.e. within a factor of 3, and optical emission varied 
within a facto r of 2. For this flare, it co uld be shown that TeV and X-ray luminosities seem to be 



connected as dCostamante et al.. 2008) 
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The l ow-state of Pks2 155-304 was monitored in a multiwavelength campaign in 2008 lAharonian et al 
d2009h . 'our wavelength bands were observed simultaneously, namely at optical (ATOM), X-ray 
(RXTE & Swift), 100 MeV to 100 GeV (Fermi-LAT) and 

> 200 GeV (H.E.S.S.) energies. While a strong correlation between optical and VHE emission 
was found, X-ray energies do not correlate with the highest energies. This behavior in the ground 
state is fundamentally different from the observed correlation between X-ray and TeV in the flar- 
ing state. The spectral index behavior observed by Fermi is around E~ 2±0 ,5 , where the deviation 
+0.5 corresponds to the temporal change in the spectral index with time. 



4. Neutrinos from Pks2155-304 

The neutrino spectrum from proton-photon interactions in blazars basically follows the primary 
proton spectrum above a brea k energy E'l? reak , where the latter is given by the kinematic condition 



for producing a A-resonance (Becker, 2008, e.g.): 



pbreak = ^ < ' K m 

(1+z) 2 80 -E 7 K ) 

Here, both the break energy E^ eak and the photon energy of the target photons, E 7 , are given in 
the observer's frame at Earth. In the equation, it was considered that the neutrino carries 1 /20th 
of the proton's energy, E v = E p /2Q. We can now express Equ. (01 numerically as 

£ break = ^ GeV . (if . (_^2_V ' . ( _L±£ \' (g) 

\10/ \MeVj \1.U6) 

Below that break energy, the neutrino spectrum is about one power flatter than the proton spec- 
trum. 

Given the observation of the low-state of Pks2155-304 in 2008 Aharonian et al. (2009), we can 



calculate the optical depth for photo-hadronic interactions. Here, we investigate if the observed 
spectral energy distribution (SED) of photons is a good target field. Using that the frequency of 
the R-band photons is v = 4.6 • 10 5 GHz, we can calculate the luminosity of the source at a given 
frequency v from the observed flux, F y , and from the distance of Pks2155-304, di, is given by 
the blazar's redshifjE 

^ = 3 ■ 1045 ' ( in io F \ , 2 ) er 8/ s • W 
\10 u 'erg/s/cm / 

Here, the reference value corresponds to the flux in the R-band, v = 4.6 - 10 GHz. Given Equ. (0, 
we can derive the optical depth for proton-photon interactions during the low-state observation 
of Pks2155-304 for the multiwavelength data. 



3 we use a ACDM cosmology with (£2 m , Q A , h) = (0.3, 0.7, 0.7) 
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Figure Q] s hows three panels: the upper one shows the SED as measured in the low-state 



( Aharonian et al. , 20091) . We can now investigate how efficient proton-photon interactions at 



different wavelengths can be and what the neutrino spectrum is. The second panel of Fig. Q] 
shows the optical depth for proton-photon interactions. At optical wavelengths, the optical depth 
is close to unity (around r py ~ 0.7), and it decreases significantly towards higher energies. Thus, 
we conclude that proton-photon interactions are efficient with optical photons, but not with higher 
energy photons. This implies that the X-ray luminosity of the source does not need to be corre- 
lated with the optical and high-energy emission, just as it is observed in the low-state. The lowest 
panel of Fig.[T|displays the break energy of the neutrino spectrum for photohadronic interactions 
for the given wavelength. It is shown that at optical wavelengths, the break energy lies around 
~ 10 8 GeV and this value decreases towards higher energy target photon fields. This is a general 
result, since the neutrino break energy is inversely proportional to the target photon field energy. 

The temporal dependence of the photo-hadronic optical depth with R-ba nd photons, derived 
from the low-state measurements of Pks2155-304 dAharonian et al. , 2009), is shown in Fig. [2] 
The upper panel shows the observed TeV emission, which correlates to the emission in the R- 
band (middle panel). The photo-hadronic optical depth is displayed in the lower panel. The 
optical depth varies in coincidence with the R-band signal, between 0.6 < T /r/o?l < 0.8, yielding 
a reasonable efficiency for the production of VHE photons and neutrinos. 

Given an average flux in the R-band of ~ 1.1 ■ 10~ 10 erg/cm 2 /s, the mean optical depth is 
o- p y opl ~ 0.7. Considering B- and V-bands would increase the optical depth further. With these 
values close to unity, the observed flux of VHE photons can be due to proton interactions with 
optical photons. The observed correlation between R-band and TeV emission can therefore be 
interpreted as a result of photo-hadronic interactions. 

Given an optical depth of around unity, and assuming VHE photons to come from photo- 
hadronic interactions, we can calculate the resulting neutrino flux for the low-state of Pks2155- 
304 using Equ. 10} . Here, we use the Fermi measurements as a reference flux. The reason for 
not using H.E.S.S. measurements is that here, we can expect to see effects like absorption, while 
the spectrum as observed by Fermi should be unaffected. The spectral behavior as observed by 



Fermi, close to dN v /d E v oc £~ 2 , seems t o mirror the primary spectrum and can therefore be in- 
terpreted as unaffected lMeli et a l. (2008); Be dnarz and Ostrowskil dl998l) . The resulting neutrino 
flux is given as the black, solid line in Fig. [3] Also shown is the predicted 1-year sensitivity for 
KM3NeT, a cubic kilometer sized neutrino detector planned to be instrumented in the Mediter- 
ranean Sea. While the energy output in neutrinos is in general high enough to be observed by 
large volume neutrin o telescopes, the ch allenge will be to have good sensitivities at energies as 
high as 10 8 GeV. In jDumm et al. , 2009b . it is reported that also IceCube is sensitive to sources 
in the southern hemisphere, if only considering the highest energy events and thereby rejecting 
the large background of atmospheric muons. IceCube's sensitivity of the 40-string configuration 
is around 

^ =5 .10-«GeV/^(^f (10) 

in the energy range of ~ 10 6 GeV -10 9 GeV. The sensitivity at these high energies is very well- 
suited for the observation of Pks2155-304. In order to observe the low-level state, the sensitivity 
needs to be increased by a factor of 10, which seems possible considering that Equ. dTOb refers 
to 330 days of the IceCube 40 string configuration. The final IceCube detector with 86 strings, 
expected to run around 10 years or more, is expected to yield significantly improved sensitivities. 
This implies that such a flux from Pks2155-304 should be detectable within the first years of 
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Figur e 1: (a) Spectral energy distribution of photons from Pks2155-304 in the low state - data from jAharonian et all 
2009). Open squares repres ent the flux at T eV energies without absorption due to the extragalactic background light, 
using the model by Kneiske and Dole ( 2008) for deconvolution. (b) Photohadronic optical depth using the SED as the 
interacting photon field, (c) Break energy of the neutrino spectrum for proton interactions with photons at a given energy. 
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Figure 2: (a) TeV lightc urve observed by H.E.S.S. iAharonian et alll20O9l ). (b) R-band lightcurve observed by ATOM 
lAharonian et al., 2009). Predicted variation of the photohadronic optical depth with time. 



IceCube. 
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4.1. Flux variability in high-energy photons and neutrinos 

In the previous section, we investigated the possible neutrino flux from Pks2155-304 in the 
low-state. This component should be present permanently. The extreme variability of this source 
at the highest energies indicates that, in addition to this permanent component, one can expect an 
enhanced contribution during the time of flaring states. We investigate this fac t at the example 



of the giant flare of Pks2155-304 observed by the H.E.S.S. experiment in 2006 lAharonian et al 



|2007l). Here, the VHE flux of Pks2155-304 showed increased activity by a factor of ~ 10 for 
about one hour. If we assume that the VHE emission is due to photo-hadronic interactions, we 
can calculate the corresponding neutrino flux according to Equ. H] The result is shown as the 
dotted lines in Fig. [3] The neutrino emission, just as the VHE photon component, is higher by 
approximately a factor of ten compared to the low-state. 

Even in this flaring state, the increased intensity can arise from an enhanced optical emission. 
On the other hand, X-ray lumi nosities are obse r ved to be coupled to the VHE luminosity with a 



correlation of Lvhe k L\ m iCostamante et alJ (120081) . The question is how to explain this very 



strong dependence of the VHE emission with the X-ray signal. O ne approach is the production 



of VH E radiation in several loops of interactions as discussed in Kelle rmann and Paulinv- Toth 
(1969) l or Inverse Compton scattering. The VHE luminosity is assumed to be produced by a 
series of n loops of interactions between the primary electrons and the secondary synchrotron 
photons, 

n 

Lvhe = L sync i, ■ ^ , (11) 
(=i 

Here, L svnc i, is the synchrotron radiation field and £ = LvHE,o/L smc h- If the initially produced 
VHE radiation Lvhe,o is larger than the synchrotron field, the series diverges for infinite loops 
n — > oo. However, if the VHE radiation is cutoff at some point when going into the Klein- 
Nishina limit, we have a maximum number n which can be small. If this happens at n = 3, 
we can reproduce the cubic correlation observed between VHE and synchrotron emission for 
the case of n = 3 and £ > 1. A similar effect could be present if the VHE radiation was not 
due to Inverse Compton scattering but due to n° decays. A detailed study of hadronic loops and 
combinations of leptonic and hadronic loops is in preparation. 



4.2. Flaring and permanent neutrino emission states 

As discussed above, an enhanced flux of photons can arise from photo-hadronic interactions, 
leading to the coincident production and emission of neutrinos. In the giant flare, the flux increase 
was about an order of magnitude for the duration of ~ 60 min. Note, however, that this enhanced 
flux is only present in the one hour when the source is flaring at very high photon energies. We 
can perform a simple exercise at what point flares can be more significant in neutrino detectors 
as KM3NeT and IceCube by considering the significance of a detection, <x ~ N S i g / ^Nbg- Here, 
N S j g is the number of signal events, while Nbg is the number of background events. The number 
of signal events in a time interval At scales as N S i g cc A v ■ At , with A v as the intensity of the 
neutrino signal, which we assume to be variable. The number of background events simply 
proportional to Nbg x At , as the background does not change with time, apart from statistical 
fluctuations. This results in a significance proportional to o~ oc Ay VAf.If we compare the flaring 
state to the permanent state, we have AjJ aie » x-Ay 61 " 1 , with A^ 316 and A^ 6 ™ as the intensities of the 
flaring and permanent neutrino states, connected by an intensity factor of x « 10. The duration 
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of the flare is Af flare « 3600 s , and Af perm is the observation time for a permanent flux. Now, we 
can compare the significances of the flaring and the permanent state, o-flare and cr perm : 



Oflare 4^ / A^e _ / Af flare 

^perm A(f™ ' V AfP e ™ * ' V AfP e ™ 1 ' 



io~ 3 



1 10/ \: 



3600 s/ \ lyr / 



This implies that the significance of a flaring state is a factor of 10~ 3 lower compared to the 
observation of a permanent flux for one year, considering that the flare's intensity is a factor of 
10 higher than the permanent emission. This leaves room for different conclusions on whether a 
flaring state or the permanent neutrino emission would be detected first: 

• Consider that Pks2155-304 flares several times a year. If we assume a constant flare activ- 
ity of x = 10, the flaring state needs to be present more than 3% of the time in one year - 
in that case, the flaring state can be observed before the permanent state. If the flaring rate 
is less, the permanent flux is likely to be observed earlier. 

• Single flares may be difficult to observed given the numbers above, unless there is no 
permanent emission, i.e. Af flare = Af perm . The optical depth we present is close to unity, 
and only a first-order approximation. If this optical depth drops far below unity most of the 
year, no permanent emission would be present and flaring states are likely to be observed 
first. In that case, we have o-flare = x • o" P erm, and selecting flaring states enhances the 
significance by a factor of x. 



5. Results and conclusions 

Very high-energy emission is observed from a growing number of blazars. One possible emis- 
sion scenario for the VHE component is photon production in photo-hadronic interactions. In 
this paper, we investigate this possibility by estimating the optical depth for proton-photon in- 
teractions in the first strong shock of Pks2155-304. Using optical photons, the optical depth 
in the low-state of Pks2155-304 is close to unity, when assuming interactions with the optical 
photon field. The observed correlation between optical and VHE wavelengths can therefore be 
interpreted as a result of photo-hadronic interactions - an increased optical flux leads to an en- 
hancement of the proton-photon optical depth and therefore to an increase of the VHE photon 
and neutrino fluxes. At the same time, no correlation between VHE and X-ray emission is nec- 
essary, since the optical depth for photohadronic interactions at X-ray energies is too small to 
lead to significant interactions. With IceCube extending the analysis of muon neutrinos to the 
southern hemisphere, the prospects of detecting the flux from Pks2155-304 are good within the 
first years of operation with IceCube. While the sensitivity of Km3NeT will in general be better, 
IceCube will be more sensitive to the important energy range, e.g. 10 8 GeV < E v < 10 10 GeV 
and both instruments would work complementary. 

Flaring states provide the possibility of enhanced neutrino emission, if the VHE radiation is 
due to 7r°-decays. These flaring states may in some cases be the only times of neutrino emission, 
if the low-state provides too low optical depth for photo-hadronic interactions. While Pks2155- 
304 provides reasonable optical depth for neutrino production also in the low-state, this may 
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Figure 3: Expected neutrino flux from photo-hadronic interactions in Pks2155-304. The solid line is the predicted 
average neutrino flux in the low-state as observed in 2008, using an E~ 2 primary sp ectrum. The dott ed line shows the 
expected emission level in the high-state of Pks2155-304 during the flare observed in Aharonian e t alj J2007h . 



not be the case for other blazars. The stacking of neutrino data during observed flaring states 
gives the unique opportunity to reduce the background of atmospheric neutrinos in VHE neu- 
trino telescopes like IceCube and KM3NeT: Even a few neutrino events from such flares will be 
significant. 
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